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Abstract 



We construct sequences of "field theoretical" lattice topological charge den- 
sity operators which formally approach geometrical definitions in 2-d CP^^^ 
models and 4-d SU (N) Yang Mills theories. The analysis of these sequences of 
operators suggests a new way of looking at the geometrical method, showing 
that geometrical charges can be interpreted as limits of sequences of field the- 
oretical (analytical) operators. In perturbation theory renormalization effects 
formally tend to vanish along such sequences. But, since the perturbative 
expansion is asymptotic, this does not necessarily lead to well behaved geo- 
metrical limits. It indeed leaves open the possibility that non-perturbative 
renormalizations survive. 

PACS numbers: ll.lO.-z, ll.15.Ha, 75.10 Hk 
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Typeset using REVT^ 



I. INTRODUCTION 



The investigation of the topological properties of 4-d SU (iV) gauge theories requires non- 
perturbative calculations. Lattice techniques represent our best source of non-perturbative 
calculations. However investigating the topological properties of QCD on the lattice is a 
non-trivial task. Topology on the lattice is strictly trivial, and one relies on the fact that 
the physical topological properties should be recovered in the continuum limit. 

In the so-called field theoretical method one constructs a local analytical function 
of the lattice fields which has the topological charge density q{x) as its classical con- 
tinuum limit: 

qL{x)^a''q{x)+0{a''+^). (1) 

qL{x) is not unique, indeed infinitely many choices differing for higher order 0{a'^^^) terms 
can be conceived. In this quite straightforward approach the major drawback is that related 
physical quantities, such as the topological susceptibility, can only be obtained after remov- 
ing cut-off dependent lattice artifacts. The classical continuum limit must be in general 
corrected by including a renormalization function 0, 

.a'^Z(<72)g(x) + 0(a^+i), (2) 

where Z{g1) is a finite function of the bare coupling going to one in the limit 0. 
The relation of the zero-momentum correlation of two qiix) operators, 

XL = Y.^qL{x)qLm. (3) 

X 

with the topological susceptibility xt is further complicated by an unphysical background 
term, 

XL{gl)=a''Z{gl)\, + M{gl), (4) 

which eventually becomes dominant in the continuum limit. So in order to extract Xt from 
this relation one needs to evaluate 2'((7q) and M((7q), which is a hard task. 

In order to overcome the problems caused by the above cut-off effects, other methods 
have been proposed. The so-called geometrical method meets the demands that the 
topological charge on the lattice have the classical correct continuum limit and be an in- 
teger for every lattice configuration in a finite volume with periodic boundary conditions. 
In 4-d SU{N) gauge theories this can be achieved by performing an interpolation of the 
lattice field, from which the principal fibre bundle is reconstructed. Since on the lattice 
each configuration can be continuously deformed into any other, integer valued geometrical 
definitions cannot have an analytical functional dependence on the lattice field. Due to their 
global topological stability, geometrical definitions should not be affected by (perturbative) 
renormalizations. On the other hand, as a drawback of their non-local nature (leading to 
non-analiticity), they may be plagued by topological defects on the scale of the lattice spac- 
ing, dislocations whose non-physical contribution may either survive in the continuum 
limit (as in the SU{2) gauge theory with Wilson action 0), or push the scaling region for 
the topological susceptibility to large (3 values. 
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Field theoretical and geometrical definitions may appear as two quite different approaches 
to the problem of studying physical topological properties. An interesting connection be- 
tween them can be conceived by exploiting the following consideration. Noting that geomet- 
rical definitions Qg can be written as a sum of local non- analytical terms, Qg = J^xQgi^)^ 
one can construct sequences of field theoretical (analytical) definitions qi^\x) formally ap- 
proaching qg{x) a.s k oo. An example of such sequences has been already given for 2-d 
CP^-i models §. 

The construction and analysis of these sequences lead to a new interpretation of the 
geometrical charges. They can be considered as limits of appropriate sequences of analytical 
operators, whose perturbative renormalizations formally tend to vanish along the sequences. 
The asymptotic nature of the perturbative expansion, which emerges from explicit calcula- 
tions of the renormalizations, does not necessarily lead to well behaved geometrical limits. 
It indeed leaves open the possibility of a residual non-perturbative renormalization, which 
would manifest itself in the contribution of lattice defects (the so-called dislocations) to the 
geometrical charge. This interpretation may clarify the meaning of the geometrical charge 
from a field theoretical point of view, and account for possible discrepancies with other 
approaches. Moreover, with increasing k, one may hope that q^^^ and their correlations 
enjoy a progressive suppression of the renormalization effects, while for not too large k their 
sensitivity to lattice topological defects is reduced. 

Our study is mainly done in the context of 2-d CP^^^ models [|,^, which are very 
useful theoretical laboratories for testing non-perturbative numerical methods conceived to 
study topological properties. We present an analysis based on large-A^ and perturbative 
calculations, and numerical simulations. We then sketch an extension to QCD, where we 
shall construct a sequence of analytical topological charge density operators approaching 
Liischer's geometrical definition |^. 



II. SEQUENCES OF TOPOLOGICAL CHARGE DENSITY OPERATORS IN 
TWO-DIMENSIONAL CP^-^ MODELS 



A. Topology in two-dimensional CP models 



Two-dimensional CP^ ^ models are defined by the action 

S=— j (fxD^zD^z (5) 

where 2 is a A^-component complex scalar field subject to the constraint zz = 1, and the 
covariant derivative = d^+iAfj_ is defined in terms of the composite field = izd^j^z. Like 
QCD, they are asymptotically free and present non-trivial topological structures (instantons, 
anomalies, 9 vacua). A topological charge density operator q{x) can be defined, 

q{x) = ^e^^d^A^, (6) 
with the related topological susceptibility 
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Xt= d^x{qix)qiO)). (7) 



A pleasant feature of these models is the possibility of performing a expansion. A 
large- analysis of the topological susceptibility shows that Xt = 0{1/N). One indeed 
finds 

.2 1 0.3801\ ^ / 1 \ 

where is the second-moment correlation length. 

Most lattice studies of 2-d CP^~^ models have been performed using the following 
actions: 

SL = -N(3Y.\z^+^Zn\\ (9) 

where (3 = 1/(2/), Zfi IS cL A^-component complex vector, constrained by the condition 
ZfiZyi = 1, and 



'S'jf^ — —N(3 ^ (^Zn+fj,Zn^n,fi + ZnZn+fi^n,fi ~ 2^ , (lO) 



where, beside the field z, a U{1) gauge field A„_^ has been introduced, satisfying Xn^^Xn,^ = 1- 
The latter lattice formulation, S\^\ turns out to be particularly convenient for a large- A^ 
expansion (see Ref. for a review). 

The original geometrical construction for the topological charge proposed by Berg and 
Liischer is 

Qg ~ Qny 

n 

Qn = TT-Ini [lnTr(P„+^+i,P„+^P„) + In Tr(P„+,,P„+^+i,P„)] , 7^ z/, (11) 

where Pn = Zn® Zn and the imaginary part of the logarithm is to be taken in (— vr, vr). For 
the lattice formulation ( |lOl) an alternative geometrical definition can be given in terms of 
the "gauge" field A„^^. Introducing the plaquette operator 

WA,n = X.^ l 2-^„_|_2,iA„ 25 (12) 

one defines Qg^\ by 

n 

ux,n = exp{i27iqx,n), (13) 



1 1 

2' 2 

volume limit of Qx „ in the form 



where qx^n ^ (— |, |). In view of large- A^ and perturbative calculations, we write the infinite 



Q\,n fiui^^n,/! ~l~ ^n+fi,i/ ^n+u,fi ^n,v)i ('^^) 
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where ^n,/^ is the phase of the field A„,^, i.e. A„,^ = e*^" ^'. 

On a finite volume and for periodic boundary conditions g„ and q\^n generate integer val- 
ues of the total topological charge for each configuration. They are not analytical functions 
of the lattice fields z and A, and fail to be defined on a zero measure set of "exceptional" 
configurations. The main feature that makes these functions different from the ones which 
we term "analytical" is the absence of single-valuedness. 

Many Monte Carlo studies of the topological properties of 2-d CP^^^ models have been 
presented in the literature. A wide range of values of has been considered, both small and 
large in order to test large- calculations. The present state of art is briefiy summarized in 
the following. 

For N = 2, that is for the 0(3) non-linear a model, recent simulations using the so-called 
classical perfect action and for a relatively large range of values of f3 seem to favour 

what suggested by semiclassical arguments, that is that Xt would not be a physical quantity 
for this model, in that a non-removable ultraviolet divergence affects the instanton size 
distribution. On the lattice this would manifest itself in the fact that lattice estimators of Xt 
do not properly scale approaching the continuum limit. Earlier claims of scaling [|I4| were 
probably originated by the relatively small range of values of /9 considered, and therefore 
scaling violations could easily be hidden behind the statistical errors. 

For N > 2 Xt should be a physical quantity. At A^ = 4 the apparent failure of the 
geometrical method is probably caused by a late approach to scaling of the geometrical 
definition (|Tl]) when using standard lattice actions such as (||) and (p!OD [p!5H T7|j7[] . Evidences 
of scaling and consistent results have instead been obtained by employing the field theoretical 
and cooling methods using action ( |T0| ) and its Symanzik improvement [0], and a geometrical 



method in the context of the classical perfect action approach |18 



At larger- A^, N ^ 10, the geometrical estimator ( [TT| ) shows scaling already at reasonable 
values of the correlation length using both action (|10]) and its Symanzik improvement []7|,p!9|. 
The available Monte Carlo data at A^ ^ 10 show that the dimensionless quantity Xt^c ^P" 
proaches, although slowly, the large- A^ asymptotic behavior (|^), and quantitative agreement 
(within a few per cent of statistical errors) is found at A^ = 21. 

In the following we shall construct sequences of operators approaching the geometrical 



definitions (O) and (13), and study their main features. 



B. Sequences of operators approaching geometrical definitions 

To begin with, we consider a sequence of analytical operators q^^n approaching the ge- 
ometrical definition g^.n expressed in terms of A„.^ fields, and within lattice formulations 
containing explicitly the gauge field, like (p^Ol). By taking appropriate combinations of the 
plaquette operator ux,n (cfr. Eq. ([T^)), one can define a sequence of local operators 
which differ from g^.n by higher and higher order terms |^ 

1 (-1)' + ^ ( 2k\ 2 , „ 

= ?A,. - p^(2-)^^?:^ + 0{qi::\ (15) 



5 



Similarly, one may define 



Ui, 



~ iTrfP P ~ ~P)V ^ ' 



and the quantities qj^n (j = 1, 2) by 



u 



exp(i27rgj,„) (17) 



(gj,„ G (— |, |)), so that one can write 

Qn = qi,n + q2,n- (18) 



A sequence of operators approaching the geometrical definition (|TT]) can then easily con- 
structed by 

Jk)_ (fc) , (k) 
Hn — Hl,n + H2,n^ 

,t) 1 * (-!)'+! I 2k\ 2 , 



(?) 

a.„-p^(2^Ar' + ofeS."')- (19) 



Strictly speaking these functions are not analytical everywhere (they are not polynomials in 
the fields) , but preserve the property of single- valuedness. 

Once the above definitions are given, one can consider sequences of lattice topological 
susceptibilities, either in terms of A or 2; fields: 

= E(?f ^S'^) = {q^''\p)^''\'P))- (20) 

Under some general assumptions, which essentially amount to assume the applicability of an 
OPE to the correlation qL{x)qL{y) when x ^ y, the relation between x^''^ and the continuum 
topological susceptibility Xt is 

X^'^\P)=a'Z^'\PYx, + M^''\(3), (21) 

where M*^'^^ can be written in terms of the identity / and the trace of the energy-momentum 



tensor operator T = yttD^zD^z, plus higher-order terms in a 120 



2/2 -^M^-^/i 



M(^)(/?) = P^'^KP) (I) + a^A^''\P) (T) + O(a^). (22) 

The expression of the background term M in terms of the identity and of the trace of 
the energy-momentum tensor (which are the only RG invariant operators of lower or equal 
dimensions sharing the same quantum numbers of Xt) may be obtained by comparing the 
OPE's of qL{x)qL{y) and q{x)q{y), and by taking their difference. 
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C. Large-AT and perturbative analysis 



(k) 

The sequence q\ ^ is particularly suitable for analytical calculations, which will show the 
main features of such constructions. The reasons that make things easier are the following. 
In both and perturbative expansion the variable 6n^^ (defined by A„^^ = e*^"-'') is used 
as a fundamental field, whose propagator can be easily derived [|1^]. In the infinite volume 
limit the expression of the geometrical charge density q\^n is linear in terms of 6n^^ (cfr. 
Eq. (0)), therefore one can easily obtain also the propagator of gA,n- The use of the g^.n 
propagator simplifies considerably the study of the large- iV behavior of the renormalization 
effects in the sequence g^'^^, which can be written as polynomials in g^.n (cfr. Eq.(|T3p). 



Indeed the leading non-trivial orders of the corresponding Z^\/3), P'^\f3), and A'^^fS) can 
be obtained by evaluating tadpole-like diagrams, whose lines are the gA,n propagators, and 
whose vertices are the coefficients of the powers of gA,n in Eq. (|T5|) . This occurs also for the 
evaluation of the leading non-trivial order in standard perturbation theory. 

In the large- limit one can unambiguosly identify the terms associated to Z^''\P), 
P^^\(3), and A^'^^P). It has been explicitly demonstrated that Qg^x and the corresponding 
susceptibility x\ the infinite volume limit are not subject to renormalizations in the 
large- iV limit |21], i.e. 



(23) 



in the continuum limit (the same has been shown for the definition ([TID too). Using Eqs. 
and (0), one can deduce the following relationship valid in the lowest non-trivial order of 
expansion 



(fc) 



where ak 



n 

{2k+l)\- 



2(2A; + l)a,(g^fo)l+«lE(^Ar^Ar) 



(24) 



-(27r 



An asymptotic expansion in powers of the large- mass |Tl|, mg, allows one to separate 
the mixing with the identity from the term of dimension two (i.e. proportional to a^) 
in Eq. (0), which contains the physical signal. This is achieved by expanding the gA,n 
propagator in powers of mo: 



(gA(p)gA(-p)) = 1^0 + mlD^ + O(m^). 



(25) 



Using results of Ref. [rT|, the calculation of the functions and is quite straightforward. 
Then by comparing Eq. (p^) with Eqs. (pT]) and (0), one can write down explicit analytical 

~ In 



Z'^\(3), and A^^\(3) in the lowest non-trivial order of 1/A^ [p2 



expressions for p[''\[3), 
the following we shall discuss some of their properties. We shall not report their complete 
expressions because they would not be very illuminating. 

The zero dimensional mixing with the identity, which dominates the Monte Carlo signal 
in the continuum limit, turns out to be 



1 



(26) 
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The renormalization functions corresponding to the terms of dimension two turns out to be 



and 



a['\p) = o 



f 1 



(27) 



(28) 



In order to derive the lowest non-trivial order contribution to A^^\l3), we have used the fact 



that in the expansion (T) = 0(1), indeed one finds [|TT 



1 

3^ 



O 



(29) 



^r(/5), 



and A^^\l3) reduce 



In the (3 ^ oo limit the large- expressions of P^\(3), Z^^^ 
to formulas that could have been obtained by standard weak-coupling perturbation theory. 
The correspondence between large- iV expansion and standard perturbation theory is obvious 
for Z^'^^jS) and P^'^\l3): it suffices to recognize that for rriQ = and /5 — > cxd the large N 
^-propagator reduces to the corresponding perturbative propagator. This correspondence 
has been explicitly verified also for A^''\P). One obtains 



Z('=)(/3) = l-k\ 



and for large k 



I3N ^ 



+ 



1 



I3N 



k+l 



(30) 



p('=)(/3) 



1 



(3N 



1 



2k+l 



(2A; + 1)!3 \(3N ^ 

> 2fc+l 

A^''\(3) ~ (2A;)! ' 



O 



BN ^ 



2A:+3' 



2k+l 



{2k)\ 



13 N ^ 



(31) 
(32) 



(fe) 



Equations ( pO| - p^ ) show the mechanism of systematic improvement of the local operators 
. As — oo, — > q\ and the lowest-order renormalizations are proportional to higher 
and higher powers of 1/(3. However, the coefficients of the leading non-trivial term grow so 
fast with k that the convergence to zero of Z('=)(/5)-l, P^^\(3) and A^^\(3) cannot be uniform 
for /3 — >■ cx). This fact leaves open the possibility that for fixed (3 some non-perturbative 
renormalization effects may eventually survive as — > oo, i.e., as the sequence approaches 
the geometrical definition. 

Inspired by this phenomenon, we suggest the following general picture. The geometrical 
charge can be interpreted as the limit of a sequence of field theoretical (analytical) operators. 
As far as (at a given (3) the renormalization effects tend to vanish along the sequence, the 
geometrical object provides a well-behaved lattice estimator of the topological charge. If 
on the contrary the renormalization effects do not disappear for k ^ oo, some pathology 
should arise in the geometrical method, such as the contribution of short-distance topological 
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defects. The asymptotic nature of the perturbative expansion, which manifests itself in 
the growing of the corresponding coefficients, leaves open the possibility of a background 
term behaving, for example, as ~ exp(— c/?), which does not get suppressed in the limit 
k ^ oo. According to the value of c, either this term is suppressed in the continuum 
limit, or it survives and spoils the expected asymptotic behavior of Xt, which should behave 
as a quantity of dimension d. The absence of non-perturbative effects at finite /5 would 
be guaranteed by a convergent perturbative expansion of the renormalizations. A good 
continuum limit of the geometrical definition would be assured by a perturbative expansion 
in which the limit k —>■ oo commutes with the continuum limit /5 — > cxd. 

The calculation of the renormalizations is more involved for the sequence expressed in 
terms of the z fields, but the conclusions are qualitatively the same. In particular, since the 
expansion of g„ in terms of perturbative fields (for example in Valent's parametrization p3| ) 
starts with a bilinear term, it is easy to see from Eq. (|15D that Z'^^\(3) — 1 and P^''\P) will 
get the leading contribution from diagrams respectively with 2k and 2k + 2 loops. Therefore 



O 



1 



2k' 



N(3 



(33) 



and 



p«(/3) = O 



2k+2 



Nj3 



(34) 



These results hold for both lattice actions (pf) and (ITO). 



D. Numerical analysis by heating method 

Since estimators involving the field A„^^ are subject to large fiuctuations in Monte Carlo 
simulations, in our numerical analysis we shall consider the sequence approaching the ge- 
ometrical definition in terms of z fields only, i.e. q^^ defined in Eq. (|19D . We present a 
numerical investigation of the corresponding renormalization effects by using the so-called 
heating method |]2^j25[| (see Ref. for an implementation of this method to 2-d CP^~^ 
models). The study of the renormalization effects by heating method relies on the possibility 
of somehow separating the various contributions in Eq. (|2l|) in off-equilibrium simulations. 
This is made possible when they are originated by lattice modes which behave differently 
under local thermalization (large vs small scale modes, Gaussian vs topological modes). As- 
suming such a distinction of lattice modes, estimates of the multiplicative renormalization 
are obtained by measuring Q^*^-* = I]n In'' ensembles of configurations constructed by 
heating an instanton-like configuration (carrying a definite topological charge Q'^^'') for the 
same number n„ of local updating steps. The plateaus showed after a few heating steps by 
data plotted as function of Uu give the desired estimates of Z'^^\ Similarly, the background 
sig nal MC^) can be estimated by measuring x^^^ ensembles of configurations constructed 
by heating the fiat configuration. Again, the plateaus of data as function of riu-, if observed, 
should provide estimates of M'^^'^ p5|,p6 . 
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We performed our simulations on a 100^ lattice. This lattice size should be sufficiently 
large in order to estimate Z and M even at large /3, in that they are expected to be short 
ranged (scaling with a), and therefore have very small finite size effects. On the other hand, 
in order to determine Xti oiis must take lattices with L ^ ^. In the following we describe 
the main results we have obtained for = 4 and = 10 using the lattice action s'f^ (cfr. 
Eq. (P). 

For /5 ~ 1 (where the correlation length is sufficiently large to expect scaling to hold |]^), 
the multiplicative renormalization Z^^^ turns out to be very close to one for small values of 
k already. For example for = 4 and at = 1.25 (which corresponds to ^ 28) we found 
^ 0.96, - 0.98, ~ 0.99, suggesting a smooth limit to one for A; ^ oo. For 
comparison we mention that the polynomial topological charge density definition considered 
in Ref. §, 

^i'^ = E e/^-Tr [P„A^P„A,P„] (35) 

(where l\^P{x) = l[Pn+^^ - Pn-f,]), has Z^p^ ~ 0.42 at this value of /?. At larger A^ Z(^) 
get closer to one, as expected from their dependence on A^. For example, at A^ = 10 and 
/3 = 1.0 (corresponding to - 17) we found Z^^^ ^ 0.985, Z^^) ~ 1.000, Z^^) ~ 1.0000, 
and ~ 0.35. 

When using the action (0), at A^ = 4 and /3 ~ 1.25 the geometrical definition (0) is 
affected by spurious contributions from short-ranged lattice structures, which turn out not 
to be negligible in order to determine Xt from Xg = ^(Q^) 0- This fact emerges clearly 
also by observing the behavior of Xg the heating of the flat configuration, in that the 
initial global geometrical charge (which is zero) appears to be soon modified by the local 
thermalization process even if the correlation length is rather large, as it occurs at /5 ~ 1.25. 
In the same heating process data of x^'^^ appear strongly correlated to those of Xg, indicating 
that the above-mentioned short-ranged configurations contribute somehow to x^''^ too (even 
for small k). In the analysis of x^''^ such spurious contributions may be interpreted as 
renormalization effects, but at these values of A^ and (3 we were not able to estimate them 
by using the heating method. It is possible that in this case the required sharp distinction 
of the lattice modes contributing to the different terms in Eq. (^Tj) does not occur, or it is 
not sufficient to make the heating method work. 

In a sense, the so-called lattice defect contributions to the geometrical definition Xg "lay 
be seen as the limit k oo oi the mixings in x^''\ which does not seem to vanish at A^ = 4 
and values of /3 corresponding to ~ 10^. Their effect probably disappears in the large- 
f3 limit. Indeed at (3 = 1.6, where the correlation length should be about one order of 
magnitude larger than at P = 1.25, we found no trace of short-ranged topological defects in 
the heating procedure. We got the following estimates of the mixing contribution M^'^'' to 
^(fc). ^{1) ^ 3 X 10-7^ m(2) < 5 X 10-^ M(3) < 10-^..., while for the lattice susceptibility 
associated with the operator (|35|) we found M^^^ ~ 7 x 10~^. P ^ 1.6, and therefore 
^ ^ 10^, may then be sufficiently large for the geometrical definition ( pi]) to be effective in 
the determination of Xt- We also quote some data at /3 = 2 where the correlation length 
should be about two order of magnitude larger than at /3 ~ 1.25, we found M^^^ ^ 2.5 x 10~^, 
M(2) < 2 X 10^1°, M(3) < 10-i\..., and M^p) ~ 2.9 x 10'^. Notice that in order to determine 
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Xt at /9 = 1.6 one would need to perform simulations on a very large lattice with L ^ ^, 
and therefore L > 10^ in order to avoid sizeable finite size effects. Of course the onset of 
scaling for Xt depends on the lattice action. In this respect actions better than (|T0|) may 
exist. A substantial improvement with respect to action (|iy) has been observed when using 
its Symankik tree- improvement |^. These results are consistent with the picture outlined 
in the previous subsection. 

As suggested by Eq. ([5^), things get improved with increasing A^. For N ^ 10 (using the 
action (|T0|)) the geometrical definition (0) seems to provide a good estimator of Xt already 
at reasonable values of the correlation length [0]. For example at A^ = 10, XgCc shows good 
scaling for .^g ~ 10 and for shorter and shorter at larger values of A^. At A^ = 10 and 
P = 1.0, the heating procedure provided the following estimates of the background signal 
in x'-'"^- M(^) ~ 0.5 X 10"^ M(2) < 0.5 x 10"^ M^^) < 0.5 x 10-^ M^^) < 10"^ and 
M(p) ~ 0.7 X 10"^ for the operator (H). The mixings in Eq. (|2l|) rapidly disappear with 
increasing k. The sequence of x^'^^ approaches Xg which in this case appears to be free of 
lattice artifact contributions. Taking into account that for A^ = 10 we have Xt^h — 0.017 [0, 
at (3 = 1.0 renormalization effects turn out to be very small and negligible for k = 1 already 
(in the evaluation of a'^Xt renormalizations for k = 1 lead to corrections of about one per 
cent). 

From the above discussion it seems that, at least when using the action (0), little 
practical improvement is achieved by the use of the operators defined by the sequences q^^y 
Indeed at low A^ (and at least for ~ 1) they appear to be sensitive to short-ranged lattice 
defects, whose contributions to renormalization effects seem as difficult to evaluate as for the 
corresponding geometrical definition (the heating method apparently fails in these cases). 
The use of standard polynomial or smeared field theoretical operators (smeared operators 
similar to those defined in Ref. and with the same features can be easily constructed for 
(jpN-i models) and cooling method p8| appears more convincing when this phenomenon 
occurs. When the geometrical definition provides a good estimator of topological activity, 
i.e. at large A^, the behavior of q^^^ is that formally predicted by perturbation theory. But in 
this case it would be more convenient to use the corresponding geometrical definition, whose 
use is further justified by our analysis. However when the renormalization effects do not 
tend to vanish (at fixed (3) along the sequence, we cannot take for granted that the physical 
predictions based on the use of the corresponding geometrical definition will be correct. 

In order to construct improved field theoretical operators one may put together both 
the idea of smearing the operators and that of defining sequences approaching geomet- 
rical definitions. This may be achieved by constructing sequences of operators in terms 
of smeared fields, instead of lagrangian fields. This should provide optimal field theoret- 
ical operators, i.e. with a more effective suppression of the renormalization effects, and 
corresponding geometrical limits probably less sensitive to unphysical short-distance lattice 
topological defects. We mention that for 4-d SU (2) gauge theory a geometrical charge de- 
fined on appropriate blocked variables has been tested p9|]. Unlike the original geometrical 
charges which turned out to be affected by dislocations, the blocked definition produced 
a topological susceptibility consistent with that obtained by alternative cooling and field 
theoretical methods. 
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III. CONSTRUCTION OF THE SEQUENCE FOR SU{N) YANG MILLS 

THEORIES 



We sketch an extension of the above study to SU {N) Yang Mills theories in four dimen- 
sions. Among the different geometrical definitions of lattice topological charge available in 
the literature, we consider that one proposed by Liischer in Ref. [Q. In the following we will 
heavily refer to this work. This geometrical definition can be written as a sum of local terms 
Qg = J2n Qn, where g„ is a gauge invariant function of the link variables. We show that it is 
possible to define a sequence of operators q^^ sharing the following properties. 

(i) For each k, q^'^ is a gauge invariant, polynomial function of the link variables of the 
single elementary hypercube c(n), possessing the appropriate (classical) continuum limit; 

(ii) The operators q^^^ tend, at least formally, to g„ for — > oo. 

qn is obtained by the lattice fields of the hypercube c(n) through a complicated interpo- 
lation procedure. An essential ingredient is the raising of link variables to fractional powers 
(this is the non-analitycal step in the definition), according to the following prescription. 
For u e SU{N), setting 

u = exp{igUaXa), (36) 

one may define the fractional power "U^, < y < 1, as 

= exp{iyguJaK) (37) 

(for convenience, we explicitly introduced the bare coupling costant g, in terms of which the 
perturbative weak coupling expansion is defined). 

We construct a sequence of gauge invariant analytical operators q^^ which differ from 
the geometrical definition g„ by 0{g'^~^^), i.e. 

qL'^ = qn + . (38) 

To this purpose, we consider a polynomial approximation U[j](i/) of u^, of degree j in u and 
u\ so that 

uy = uy]{y) + 0{g'^+'). (39) 

The polynomial functions M[j](y) can be easily constructed. We then make the appropriate 
substitutions '^[jliv) ^he expression of Liischer's topological charge density so to 
obtain q^^^ as an approximation 0{g'''^^) of g„ that contains polynomials of minimum degree 
in the link variables. 

The lattice susceptibilities x*-^^ = J2n{Qn^Q(t^) should be related to the continuum topo- 
logical susceptibility Xt by 

X^'^\P)=a'Z^'^\f3Yx, + M(''\(3), (40) 

where 

M«(/?) = p('=)(/3) (/) + a^A^''\l3) (T) + 0{a^). (41) 
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Assuming non-renormalization for the geometrical definition g„, and using Eq. (^) one can 



infer that {(3 = ^) 



zW(/5) = 1 + 0(^1), (42) 
P''\(^)=0(-^], (43) 



where 



/ = — for even k, 

2 

/ = for odd k. (44) 

An analysis of the behavior of the coefficients 2;*^^^ and p^'^^ of the leading non-trivial order 
m Z^^^ and P('=) can be more easily carried out for even k, for which all the corresponding 
diagrams are substantially tadpoles. A rough estimate of the behavior of p^''^ and z^^^ may 
be obtained by counting all the contractions to form k/2 tadpoles: 

^ {k/2)\ ^ ^ 

Here the dependence on N has been overlooked, in that all contractions have been considered 
as giving the same contribution. This is true only in the case of commuting generators, as 
in the case of U{1) gauge theory. Neverthless Eq. (|^) should give an idea of the behavior 
at large k, at least for not too large N. On the other hand, at large- contractions of non- 
sequential generators in the traces are suppressed by powers of 1/A^, so that the contractions 
contributing at N = 00 are considerably reduced, leading probably to p^''^ ~ z^''^ ~ 0(1). 

Then, as already observed for the sequences constructed within the CP^^^ models, the 
lowest-order renormalizations are proportional to higher and higher powers oil/ 13. However, 
at finite A^ the corresponding coefficients of the leading non-trivial order grow so fast with 
k that the convergence to zero of the renormalization functions cannot be uniform for j3 — >■ 
00. Again, with increasing A^ renormalization effects should be further suppressed in the 
sequence, suggesting that at least at large-A^ the geometrical charge should be free from 
dislocations. 
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